Exclusion limits from data of directional Dark Matter detectors 
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Directional detection is a promising search strategy to discover galactic Dark Matter. Taking 
advantage on the rotation of the Solar system around the Galactic center through the Dark Matter 
halo, it allows to show a direction dependence of WIMP events. Even though the goal of directional 
search is to identify a WIMP positive detection, exclusion limits are still needed for very low exposure 
with a rather large background contamination, such as the one obtained with prototype experiments. 
Data of directional detectors are composed of energy and a 3D track for all recoiling nuclei. However, 
to set robust exclusion limits, we focus on the angular part of the event distribution, arguing that 
the energy part of the background distribution is unknown. As the angular distributions of both 
background and WIMP events are known, a Bayesian approach to set exclusion limits is possible. 
In this paper, a new statistical method based on an extended likelihood is proposed, compared 
to existing ones and is shown to be optimal. Eventually, a comprehensive study of the effect of 
detector configurations on exclusion limits is presented. It includes the effect of having or not 
sense recognition, a finite angular resolution, taking into account energy threshold as well as some 
astrophysical uncertainties. 

PACS numbers: 95.35.+d, 14.80.-j 



I. INTRODUCTION 



In the context of direct detection of non-baryonic Dark 
Matter (WIMP), an alternative strategy to massive de- 
tectors [l|-y|, aiming at high background rejection and 
planning to scale up to ton-scale, is the development of 
detectors providing an unambiguous positive WIMP sig- 
nal. This can be achieved by searching for a correlation 
of the WIMP signal with the solar motion around the 
galactic center, observed as a direction dependence of 
the WIMP flux 3, coming from (£ Q = 90°, b Q = 0°) in 
galactic coordinates, which happens to be roughly in the 
direction of the Cygnus constellation. This is generally 
referred to as directional detection of Dark Matter and 
several projects of detector are being developed for this 
goal [I|5|. 

The major advantage of directional detection is the fact 
that the angular distribution of WIMP events is pointing 
toward the Cygnus constellation while the background 
one is isotropic. This opens the possibility to really iden- 
tify the WIMP signal as such. The Discovery param- 
eter in this case is a signal pointing toward Cygnus as 
emphasized in [lol since it cannot be mimicked by 
background. Recently, a statistical map-based analysis 
has been developed showing the possibility to ex- 
tract from data samples of forthcoming directional de- 
tectors, both the main direction of the incoming events, 
thus proving the galactic origin of the signal, and the 
number of WIMP events contained in the map thus con- 
straining the WIMP-nucleon cross section. 
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Even though the goal of directional experiments is to 
identify a WIMP positive detection, exclusion limit is 
still needed, for very low exposure with a rather large 
background contamination, such as the one obtained with 
prototype directional experiments. 

The directional insensitive direct detection strategy is 
based on the measurement of the recoil energy (dN / dEp) 
and a dedicated statistical method has been developed 
[l2j to optimize the exclusion limit obtained from a given 
set of data. In addition to the recoil energy, directional 
detectors also provide, for each track, the scattering an- 
gle (f2). Recently a 2D extension of the previous method 
has been proposed [H, [id]]. It is based on the double- 
differential spectrum (d 2 N / dEudVt) and allows to ac- 
count for all information given by a directional detector. 
As a matter of fact, the energy spectrum of the back- 
ground is unknown [l5l . [l6j , while its angular spectrum 
is expected to be isotropic in the galactic rest frame. We 
take advantage on this point by proposing a new method, 
taking into account only the angular part of the spectrum 
(dN/dCi) in a given recoil energy range. As no assump- 
tion on the background energy dependence is needed, ro- 
bust and conservative exclusion limits may be provided. 
We also argue that upcoming directional data may po- 
tentially be contaminated with a rather large amount of 
background events and this method is intended to cope 
with low signal to noise ratio, when identification is not 
possible This method is based on a likelihood anal- 
ysis of events that are supposed to be composed of both 
background and WIMP events. 

The paper is organized as follows. In section [TT| we 
will first introduce the directional detection framework. 
In section IHIl we will introduce the different statistical 
methods used to set WIMP-nucleon cross section limits 
: the existing one (Poisson, Maximum Gap [Hj]) and a 
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new one, referred to as the Directional Likelihood method. 
They will be compared in section IIVI for various detec- 
tor configurations such as ideal detector, with or without 
sense recognition, finite angular resolution, and threshold 
effect. It allows to evaluate the impact of experimental 
uncertainties on each method. The effect of some astro- 
physical uncertainties like the WIMP local velocity dis- 
persion and the local dark matter density po will also be 
discussed. 



II. DIRECTIONAL DETECTION: 
EXPERIMENTAL AND THEORETICAL 
FRAMEWORK 



Several directional detectors are being developed 
and/or operated : MIMAC @, DRIFT @, DM-TPC 
and NEWAGE [8]. A detailed overview of the status 
of experimental efforts devotedjo directional Dark Mat- 
ter detection is presented in [9]. Directional detection 
of Dark Matter requires track reconstruction of recoiling 
nuclei down to a few keV. This can be achieved with low 
pressure gaseous detectors [17J and several gases have 
been suggested : CF 4 , 3 He -I- C4H10 or CS 2 . Both the 
energy and the track of the recoiling nucleus need to be 
precisely measured. Ideally, recoiling tracks should be 
3D reconstructed as the required exposure to reject the 
background hypothesis is decreased by an order of magni- 
tude between the 2D readout and 3D readout [18( . Sense 
recognition of the recoil track (head- tail) is also a key 
issue for directional detection [19h21| and its impact on 
exclusion limits calculation will be discussed here. The 
study is done for an ideal, but realistic, detector which 
could be within reach in a few years. We consider a 10 
kg CF4 detector, operated at 50 mbar and allowing 3D 
reconstruction of recoiling tracks. A recoil energy range 
[Er 1 ,Er 2 ] is chosen between 5 and 50 keV. The lower 
bound of the energy range is due to the threshold ioniza- 
tion energy taking into account the quenching factor. As 
most of the WIMP events are concentrated at low recoil 
energy, an upper bound is chosen to limit the background 
contamination of the data. Indeed, for a WIMP mass of 
100 GeV.c~ 2 , 70% of the recoils are between 5 keV and 
50 keV and only 10% are above 50 keV in the case of 
an escape velocity taken as v esc = 00 and a form factor 
F(Er) approximated to one (45|. Thus, increasing the 
upper bound would lead to a potentially weaker signal 
to noise ratio. In section ITVl we will study the effect of 
having or not sense recognition, degrading the angular 
resolution as well as varying the energy threshold. 
Through this study, we focus on the simplest model for 
the Milky Way halo : the isotropic isothermal sphere in 
which the WIMP velocity follows a Maxwellian distribu- 
tion defined in the laboratory rest frame as 
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with a dispersion <r v = vo/y/2 where vo = 220 km.s -1 is 
the circular speed at large radii. Effect of uncertainties on 
vo, i.e on the local WIMP velocity dipersion, on exclusion 
limits will be discussed in sec. IIVEI We consider a de- 
tector velocity equal to the tangential component of the 
Sun's motion around the Galactic center vq = 220 ± 20 
km.s -1 , neglecting the Sun's peculiar velocity and the 
Earth's orbital velocity about the Sun. As a matter of 
fact, their contribution to the detector velocity is smaller 
than the uncertainty on v@ [13] • Using the galactic coor- 
dinates (£, b), the WIMP velocity is written in the galac- 
tic rest frame as : 

v — v (cos £ cos b x + sin I cos b y + sin b z) 

where £ points toward the galactic center, y in the direc- 
tion of the solar motion and z toward the galactic north 
pole. The WIMP-induced recoil events are then com- 
puted by generating random incident WIMP velocities 
from vf(v) and assuming an isotropic elastic scattering 
in the center of mass frame. Then, the recoil energy (En) 
is given, in the laboratory frame, by : 



E R = - * cos 2 9 R 



- m N ) 
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with m x the WIMP mass, tun the mass of the target and 
Or the recoil angle with respect to the WIMP direction. 
Thus, this Monte Carlo WIMP event generation can be 
used both for simulating experiments and for evaluating 
theoretical angular distributions. Then, in the case of the 
binned version of the Directional Likelihood method, de- 
scribed in sec. MI C 1[ the theoretical event distributions 
of WIMP events S and background events B are esti- 
mated using 10 8 Monte Carlo generated events. In the 
case of the background, the events are randomly gener- 
ated according to an isotropic distribution. However, we 
have also developed an unbinned definition of our new Di- 
rectional Likelihood method (see sec. IIII C ip . In this case, 
the S and B distributions are analytically calculated. In- 
deed, as the spherical isothermal halo model leads to a 
WIMP-induced recoil distribution which presents a spa- 
tial distribution of events with an axial symmetry along 
the y axis, we will use the cos 7 distribution of the events 
with cos 7 = y.r, f being the recoil direction in the galac- 
tic coordinates. In order to evaluate S(cos7) (see sec. 
IIII C 1[) . we used the following theoretical distribution 



dR 

dcos7 



k I exp 
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dE R (3) 
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k is a normalization factor, v m i n = ^Jm^ERjlm 2 
corresponds to the minimal WIMP velocity required 
to produce a nuclear recoil of energy Er and m r is 
the WIMP-nucleus reduced mass. Moreover, it is 
worth noticing that due to this axial symmetry, there 
is no loss of information when looking at ID angular 
distribution (cos 7) rather than 2D angular spectrum 
(I, b). The theoretical background distribution £?(cos7) 
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corresponds to a flat distribution as the background is 
taken to be isotropic. 

The expected number of WIMP events p is calculated 
by integrating the theoretical recoil energy distribution 
over the experimental energy range to Er 2 : 

with o the WIMP-nucleus elastic scattering cross section, 
po the local dark matter density and £ the exposure. In 
the following we define Co the WIMP-nucleon cross sec- 
tion directly related to a in the pure-proton approxima- 
tion, see sec. [V] 

III. WIMP-NUCLEON CROSS SECTION 
LIMITS CALCULATION METHODS 

In this section, we present the existing statistical meth- 
ods, Poisson and Maximum Gap .12], and we propose a 
new one dedicated to directional detection, referred to as 
Directional Likelihood method. 



A. Poisson method 

The simplest way to derive a cross section upper limit 
from a given experimental dataset, is to use the Poisson 
method. This is the most conservative approach because 
it assumes no knowledge on neither the background nor 
the WIMP event distribution. Each recorded event is in- 
terpreted as a WIMP one. Hence, the aim of this method 
is to define the value of <To corresponding to a number of 
signal events, n eX c, which is excluded at a given confi- 
dence level (CL) according to the Poisson distribution. 
The confidence level is then given by 1 — a(fj,), with 

N m 

a(A») = e-"E a T (5) 
* — ' m! 

and N is the total number of events contained in a 
dataset. Then, the excluded value of cr at 90% CL is de- 
duced from the value of /x eX c, which satisfies 1— a{jjL exc ) — 
0.9, using equation ((U). 

As in the Poisson method all events are considered as a 
WIMP signal, any background observed events will lead 
to overly estimated upper limits. Even though it has 
been shown [l2|, EH that, in the case of a pure WIMP 
signal, the Poisson method gives more restrictive upper 
limits than the Maximum Gap method, in this paper we 
are only interested in a dataset without WIMP events 
or with a large background contamination, which seems 
to be closer to real data from upcoming directional de- 
tectors Q. We argue that in the case of a pure WIMP 
signal, directional detection should lead to an identifica- 
tion of a positive WIMP detection, using the methods 
presented in [Io|, [H[ • 



B. ID Maximum gap method 

In order to deal with an unknown background contam- 
ination, which is the case with direction-insensitive direct 
detection {i.e. measuring only the energy of the recoil), 
two statistical methods have been proposed by S. Yellin 
[l2| : " the maximum gap" and "the optimal interval" . In 
this paper, we define a modified version of the maximum 
gap method applied to ID directional data, i.e. the an- 
gular part of the spectrum (dR/dcosj distribution, see 
sec. |TT}. Then, for a given value of do and a given dataset 
of N recorded events, we obtain a set of (N+l) gaps Xi 
which are calculated as 12]: 

/■cos7 i+ i dR 

Xi = / (<r )dcos7 (6) 

,/co S7i dcos7 

with COS70 = — 1 and cos 7jv+i = 1 which are the lower 
and upper bounds of the dN/d cos "/ distribution. The 
maximum gap x is then defined as the largest Xj ob- 
tained for a given experiment. The interest of the "max- 
imum gap" method is the fact that the distribution of 
x is independent of the theoretical distribution of the 
events. It depends only on fi and not on the shape of 
the event distribution. We refer the reader to [l2| for a 
detailed discussion on setting limits using the Maximum 
Gap method. 

We do not use full information from the datasets, i.e. 
d 2 R/dEnd cos 7, as in [13| . since the performance of the 
method strongly depends on the assumed background 
distributions. As the energy distribution of background 
events is unknown, we prefer to consider only the angu- 
lar part of the distribution for which the background is 
known to be isotropically distributed, leaving aside the 
energy information of recorded events. 

C. Directional Likelihood methods 

As stated above, we restrict the information to the 
angular distribution of the events. In this case, the back- 
ground is well understood and then, we have a theoretical 
distribution for both WIMP events and background ones 
thus allowing a Bayesian calculation of exclusion limits 
using a likelihood analysis. 

1. Definition 

The Directional likelihood method is based on a recent 
paper [l(| which aims at distinguishing a genuine 
WIMP signal in a directional dataset. This likelihood 
definition allows to recover the main direction of the 
recoils and the number of WIMP events contained in the 
recoil map in order to identify a detection of particles 
from the galactic Dark Matter halo and to evaluate its 
significance. It has been shown that this analysis tool 
gives satisfactory results on a large range of exposure 



and background contamination levels. But, for very low 
number of WIMP events and for very large background 
fractions, when this method failed at recognizing a 
WIMP detection, obviously an exclusion limit should be 
derived. This will be the case for the very first results of 
directional detectors with low exposures. 

The total number of recorded events N is considered as 
the sum of n s signal events and rib background events, i.e. 
N = n s + rib, where both n s and rib can be regarded as 
Poisson variables with means fi s and /if, respectively. The 
fact that both signal and background angular spectra are 
well known allows to derive upper limits using the Bayes' 
theorem. In the case of flat priors for both fi a and fib, 
and taking the evidence as a normalization factor, it is 
reduced to 

P(p s ,Hb\D) oc if (ji B ,fib) 

where D refers to the characteristics of the data, as the 
total number of recorded events N, their direction and 
their energy. 

However, in order to incorporate in the likelihood def- 
inition both the information contained in the measured 
angular spectrum and the fact that iV is a Poisson vari- 
able of mean fi^ = Ms + Mfcj an extended likelihood func- 
tion is used: 

#0*., /*) = ^+^V ^> x L(fi s , fi b ) (7) 

In the following, two different definitions of L(fi s ,fib) 
are proposed, in the case of binned or unbinned data. 

When using binned data, the following definition of 
L(fi s , Hb) is used : 

L(fl s ,fl b ) = IT P ( T Si ~ l X B i\ M i) 

7=1 V Ms + M6 Ms + Mb / 

where M corresponds to the observed recoil map, S 
and B refer to the signal (WIMP) and background 
theoretical distributions (see sec. HI]). The probability P 
is then evaluated at each bin with a poissonian distribu- 
tion. Here the bins are the pixels of the observed recoil 
map represented in galactic coordinates. An angular 
resolution of 15° (FWHM) is chosen [23], corresponding 
to a number of equal area pixels of N p i xc \ s — 768. 

Unbinned data are often preferred and we propose a 
unbinned version of the definition of L(fi s ,fib), defined 
as : 

N , v 

L(fi s ,fi b ) = TT — ^S(cos 74 ) + _J^B( C0S7i ) 
f=iVMs+Mf> Ms + Mfc / 

with S and B the theoretical distributions, dR/d cos 7 
(see sec. |H]), of WIMP and Background events respec- 
tively. 
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FIG. 1: Representation of «Sf (fi s , Mb) i n the case of expected 
WIMP events and 30 expected background events. The lower 
left panel presents the 2 dimensional likelihood distribution 
with the 68%, 90% and 95% CL contours. The orange star 
represents the input values of /t s and fib- The left upper 
panel and the right lower panel represents the marginalized 
distributions P(jj, s \D) and P(/j, b \D), respectively. 



2. Setting exclusion limits using P(fi s , fit,\D) 

Using a Bayesian approach, the probability density 
function of the parameter of interest fj, s can be derived by 
marginalizing P(/i s , fj,b\D) over the parameter \ib- Then, 
the value of /i cxc is obtained by solving: 

P(ri s \D) dLi s =0.9, (8) 

Jo 

Thus, the value of the excluded cross section at 90% CL 
is deduced from /x oxc using equation (jj). As an illustra- 
tion of the method, figure [I] presents the results from the 
calculation of if (/i s , /if,) in a working example with ex- 
pected WIMP event and 30 expected background events. 
The probability density functions of /x s and /Jb are de- 
rived from the marginalization of the 2 dimensional like- 
lihood over the /if, and /i s parameters respectively. The 
calculation of /i cxc is also illustrated and represented on 
the left upper panel of figure [TJ Hence, on this example 
of 30 simulated events, the method allows to exclude at 
90% CL a number of WIMPs greater than 7.8. At last, a 
strong correlation (p ~ —I) between the parameters /i s 
and /if, is observed on the 2 dimensional representation of 
if (/i s , fib) and is explained by the fact that n s = N — rib. 

3. Comparison between binned and unbinned likelihood 
definitions 

When using binned data, bias can be introduced in 
the analysis and thus, whenever possible, unbinned data 
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analysis are preferred. For each configuration and each 
result presented in this paper, both binned and the un- 
binned analyses have been done, using the proper like- 
lihood definition. No difference between them has been 
seen, due to the fact that, in our case, the binning (pix- 
els area) is thin enough and the theoretical distributions 
are smooth and much wider than the pixels. As the ef- 
fect of bining is null, using binned data do not introduce 
any bias in the calculation of P(/z s ,/Zb|D) and then in 
the estimation of the excluded cross section. Hence, only 
the results obtained with the binned likelihood defini- 
tion, referred to as the Directional Likelihood method, is 
presented hereafter. 



D. Discussion on limit setting methods 

In order to evaluate the frequency distributions of the 
excluded cross section at 90% CL, see figure [21 we used 
10,000 toy Monte Carlo experiments for each case and 
for the three statistical methods. 

From each distribution, we can derive the median value 
of the excluded WIMP-nucleon cross section <7 me( j and 
the RMS given by the Poisson, Maximum Gap and the 
Likelihood methods. As an example, on figure [2] we 
consider a case with 10 expected WIMP events and 10 
expected background events. The coverage is defined 
as the fraction of experiments leading to an excluded 
cross section above the input value. The values deduced 
from the study are : (99.99 ± 1)%, (98.76 ± 1)% and 
(93.39 ± 1)% for the Poisson, Maximum Gap and likeli- 
hood methods respectively. The fact that the coverage is 
greater than 90% for the three different methods means 
that in the case of a background contamination, the ex- 
cluded cross section is overly estimated. In the case of 
the Poisson method, this over estimation comes from the 
fact that each recorded event is considered as a WIMP 
event. Contrary to the Poisson method, the Likelihood 
and the Maximum Gap methods evaluates the agreement 
between the distribution of the observed events with the 
theoretical WIMP event distribution. The fact that the 
coverage of the Likelihood method is closer to 90% CL 
than the other one is explained by the incorporation of 
the L(/i s ,iib) term in the extended likelihood function 
which helps at discriminating WIMP events from the 
background ones contained in a given dataset. Thus, 
the Likelihood method is the most robust method in the 
case of a sizeable background contamination. Note that 
if the expected number of background events could be ac- 
curately estimated even with some uncertainties, it could 
be possible to consider a non constant prior P(/i&) in or- 
der to improve the effici ency of the Likelihood method. 
However, as discussed in [L^, as the number of expected 
background \Xb cannot be known without important un- 
certainties, it is not possible to use the Feldman- Cousins 
technique [2J] and the ones presented in [2a ] . 
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FIG. 2: Frequency distributions of the excluded cross section 
<r cxc at 90% CL for each method with 10 expected WIMP 
events and 10 expected background events. The Poisson, 
Maximum Gap and Likelihood methods are represented by 
the black solid line, the red dashed line and the blue dotted 
line respectively. The coverage for each methods are (99.99 ± 
1)%, (98.76 ± 1)% and (93.39 ± 1)% respectively. 



IV. COMPARISON OF THE EXCLUSION 
LIMIT METHODS 

As the statistical framework is being defined, we can 
compare the efficiency of the three methods. We focus 
on pure background or highly background-contaminated 
data, as it is the most realistic case for upcoming direc- 
tional detectors, especially prototype experiments with 
low exposures. In the following, we evaluate with the 
three statistical methods, the effect of detector config- 
urations on the exclusion limits, such as having or not 
sense recognition capability, having a low angular res- 
olution and varying the energy threshold. The effects 
of astrophysical uncertainties on the calculation of cross 
section limits are eventually evaluated. 
Unless otherwise stated, we consider a CF4 directional 
detector and a halo model as presented in sectionflll with 
a total time exposure of 3 years and a WIMP mass taken 
to be 100 GeV.c- 2 . 



A. Ideal detector 

First, an "ideal detector" is considered, meaning 
a detector allowing a 3D reconstruction of recoiling 
events with sense recognition and with perfect angular 
definition. The three methods are compared on two 
simulated datasets : pure background (fig. [3]) and 
background events with 5 expected WIMP events 
(fig. The first case reflects a scenario in which the 
exposure is not sufficient to access WIMP events due 
to a very low WIMP-nucleon cross section while in 
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the second case a few WIMP events are obtained as 
the considered cross section was increased. Further an 
increase of the number of WIMP events is meaningless 
in the context of directional detection, as in such case 
a discovery could be done rather than an exclusion, as 
shown in [10|. Figure [3] and H] present the excluded cross 
section (<T me <2 a s defined in sec. IIII D[) obtained with the 
three different methods as a function of the number of 
expected background events. 

In the pure background case (fig. [3]), the three methods 
start with a common value at null event. Indeed, in 
such a case, the Maximum Gap and the Likelihood 
methods are identical to the Poisson one. Then, when 
no event is observed, an excluded value of /i cxc — 2.3 
WIMP events at 90% CL is deduced from the three 
methods. The excluded cross section obtained with the 
three methods start being different as the number of 
background events increases. The Likelihood method 
allows to obtain a much more restrictive exclusion limit 
than the two others, by a factor ~ 4 at high background 
contamination (30 events). As expected, increasing the 
number of background events degrades the upper cross 
section limit (i.e. increases the value of u m ed)- This 
holds for the three methods, noticing that the Likelihood 
method is less sensitive to background contamination, as 
it takes full advantage on the knowledge of the expected 
WIMP and background angular distributions. 
In the case of data populated with a few WIMP events, 
the excluded cross section evaluated with the Likelihood 
method is only slightly less restrictive than in the former 
case. The two Poisson curves (black solid lines) are 
exactly the same on figures [3] and [4] Indeed, as the 
Poisson method is only dealing with the total number of 
events, the presence of WIMP events does not modify 
the exclusion limit. It can be noticed that below ~ 10 
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FIG. 3: The excluded cross section (a me d) as a function of 
the number of background events, in the case of a pure back- 
ground, obtained by the three different statistical methods: 
Poisson (black solid line), Maximum Gap (red dashed line) 
and Likelihood (blue dotted line). 
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FIG. 4: The excluded cross section (cr mec j) as a function of the 
number of background events, in the case of 5 expected WIMP 
events, obtained by the three different statistical methods: 
Poisson (black solid line), Maximum Gap (red dashed line) 
and Likelihood (blue dotted line). 



events, the Maximum Gap method gives worse results 
than the Poisson one. As a matter of fact, in the case 
of pure WIMP signal the Poisson method overcomes the 
Maximum Gap one, see sec. IIII Al 

As shown in [10], the directionality of the WIMP 
signal depends on the WIMP mass. Lighter is the 
WIMP, stronger is the angular anisotropy. Indeed, due 
to the finite energy range and the fact that low WIMP 
mass induces an energy distribution shifted to low 
energy, events above threshold (E^) are the one with 
the most directional feature (eq. ©). The directionality 
evolves quickly at low masses and very slowly for masses 
heavier than a hundreds of GeV.c" 2 . A study of the 
effect of the WIMP mass on the cross section limit has 
been done. Figure [5] presents on the upper panel the 
cross section limit (cr m ed) as a function of the WIMP 
mass for the three methods considering 10 expected 
background events and one expected WIMP event. 
This highlights the fact that the Likelihood method 
gives the best result on the whole WIMP mass range. 
Moreover, we can notice that the cross section limit is 
more restrictive for masses lower than ~ 30 GeV.c~ 2 
due to the threshold of the detector. Indeed, at low 
masses, the energy distribution is more peaked at low 
energy leading to a loss of expected events and then 
to an increasing cross section limit. The lower panel 
of figure [5] presents the evolution of <T me d divided by 
do (the input value to get 1 WIMP event). It allows 
to show that the cross section limit is enhanced by 
the directionality of the signal, at low WIMP masses 
(< 100 GeV.c -2 ), for both the Maximum Gap and 
the Likelihood method. Increasing the WIMP mass 
does not affect the result as there are little differences 
in angular distributions for heavy WIMPs, as stated 
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FIG. 5: Upper panel : cross section limit (<r me d) as a func- 
tion of the WIMP mass for the three methods considering 
10 expected background events and one WIMP event. Lower 
panel : cross section ratio (owd/co) as a function of the 
WIMP mass. This highlights the fact that the cross section 
limit is enhanced by the directionality of the signal, at light 
WIMP masses (< 100 GeV.c -2 ), for both the Maximum Gap 
and the Likelihood method. 



above. As expected the Poisson method is not sensitive 
to the WIMP mass, as it handles only the total num- 
ber of events recorded and not their angular distribution. 

A key feature of the proposed Likelihood method is the 
fact that, contrary to the Maximum Gap method [Hj], 
it can deal with very large number of events. Hence, 
we study hereafter the evolution of a me d with the expo- 
sure, for a fixed WIMP mass m x — 100 GeV.c -2 . Figure 
[6] presents the excluded cross section obtained with the 
Poisson and Likelihood methods as a function of the ex- 
posure. Curves are labeled according to the expected 
value of A, defined as the WIMP event fraction in the 
dataset (A = /j s / (/x s + /it)). To illustrate, we have chosen 
two different background rates (2, 20) evt/kg/year and an 
input cross section do = 10 -5 pb giving an expected rate 
of WIMP events equals to 0.2 evts/kg/year. The Poisson 
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FIG. 6: The excluded cross section obtained with the Poisson 
and Likelihood methods as a function of the exposure. A is 
defined as the WIMP fraction in the events : A = /Lt s /(// s +//(,). 
For this study the WIMP mass is fixed m x = 100 GeV.c" 2 
and two different background rates are considered : (2, 20) 
evt/kg/year. The input cross section is a — 10 -5 pb giving 
an expected rate of WIMP events equal to 0.2 evts/kg/year. 



results (in black circles) are also represented in order to 
appreciate the efficiency of the Likelihood method (blue 
triangles). 

With the Likelihood method, the cross section upper 
limit is still getting more restrictive, when increasing the 
exposure, especially in the case of highly background- 
contaminated data. On the contrary, if the background 
fraction is too high (A < 0.01), the cross section limit 
given by the Poisson method remains flat as a function of 
the exposure. This is the proof that the Poisson method 
will not be able to deal with data contaminated by a large 
number of background events, as expected from upcom- 
ing directional prototype detectors. Then, the difference 
between the Poisson and the Likelihood method is about 
an order of magnitude on this range of exposure. But 
the most interesting point is the fact that only the Like- 
lihood method will allow to set better exclusion limits 
with increasing exposure, even with critically low value 
of A. As expected, when increasing the rejection power of 
the detector, i.e. increasing the expected value of A, the 
two statistical methods (filled triangles and filled circles) 
are getting more efficient. 



B. Detector without sense recognition 

Even though several progresses have been done p~9i — 
H3, sense recognition (so-called "Head- Tail" (HT) ef- 
fect) remains a key and challenging experimental issue 
for directional detection of Dark Matter. In particu- 
lar, it should still be demonstrated that sense recogni- 
tion may be achieved at low recoil energy where most 
WIMP events reside. Hence, in the following we inves- 
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FIG. 7: The median upper limit cross section obtained by 
the three different statistical methods as a function of the 
number of background events in the case of pure background. 
The results without sense recognition are presented with filled 
markers while empty markers refer to a detector with sense 
recognition. 



tigate the effect on exclusion limits of using a detector 
with no sense recognition. In such case, a recoil coming 
from (cos7,</>) cannot be distinguished from a recoil com- 
ing from (— cos 7,0 + 71-). Then, we need to use axial data 
characterized by a direction (| COS7I) and the directional 
event distribution becomes: 

dN _ dN (cos + diV(-cos7) ^ 
rf|cos7| dcos7 dcos7 

We study the effect of no sense recognition in the case 
of pure background data and the result is presented on 
figure [7] To ease comparison, results with sense recog- 
nition are recalled (from figure [3]). Still, the Likelihood 
method overcomes the two others and it appears that the 
absence of sense recognition only mildly alter the result 
(a factor of three at high background contamination). 
As expected, the Poisson method is insensitive to the 
absence/presence of sense recognition. Surprisingly, the 
Maximum Gap method becomes less competitive than 
the Poisson method in the case of no sense recognition. 
This is due to the fact that in this case, the expected 
WIMP distribution becomes less anisotropic and then it 
is getting closer to the expected background event dis- 
tribution (isotropic). Indeed, as it has been shown in 
[HI, [l3j], in the case where the difference between the 
background and the WIMP event distribution is very 
small, the Poisson method sets more competitive limits 
than the Maximum Gap one. 

Taken at face value, this result suggests that sense recog- 
nition may not be so important for directional detection 
when trying to set exclusion limits. The difference be- 
tween 100% sense recognition on the whole recoil energy 
range, which is obviously unrealistic, and no sense recog- 
nition is only minor (less than a factor of three at high 



background contamination). The worst case is indeed a 
partial sense recognition strongly depending on the re- 
coil energy. In this case, we suggest not to consider this 
information to set exclusion limits. However, sense recog- 
nition remains a key issue worth investigating, for WIMP 
discovery which is the ultimate goal of directional detec- 
tion. 



C. Detector with finite angular resolution 

Previous studies have been done in the case of a de- 
tector with a perfect angular resolution. However, real- 
istic data of upcoming directional detectors should suffer 
from finite angular resolution. This is an intrinsic limita- 
tion of this detection strategy. Even if simulations show 
that straight line tracks may be 3D reconstructed with 
a rather small angular dispersion [26|, realistic tracks in 
low pressure gaseous detectors would encounter a rather 
large "straggling" effect (angular dispersion). The lower 
is the recoil energy, the larger is the angular straggling. 
Hence, in the following we investigate the effect on ex- 
clusion limits of using a detector with a finite (realistic) 
angular resolution. Having a finite angular resolution 
means that a recoil initially coming from the direction 
f(8,(p) is reconstructed as a recoil f' (6', (/)'). Then, the 
angular deviation between the initial recoil and the re- 
constructed one is given by: 

e = cos _1 (f . f) (10) 

The f' direction is characterized in the f frame by the 
angles and $ which are, for each event, randomly gen- 
erated using the following probability density function 

/(e,$) : 

/(e,$)ocsmeexp^-^j i $ i 2tt W 

In order to compute the theoretical angular distribution 
considering this finite angular resolution, we have done 
the convolution product of the initial angular spectra 
dN/ dQ, with a smoothing gaussian function of width o& 

Figure[5]presents the cross section limit as a function of 
the angular resolution hq in the case of 10 expected back- 
ground events and one expected WIMP event. It can be 
noticed that the Maximum Gap and the likelihood meth- 
ods are only slightly dependent on the angular resolution 
of the detector. The deviation for the two methods is of 
the order of 30 % from <jq — 0° to uq — 45°. Hence, as 
far as exclusion limits are concerned, the effect of angu- 
lar resolution is relatively small but has to be correctly 
evaluated to set coherent upper limits. This outlines the 
need for detector commissioning, e.g. by using a neutron 
field [H. 

As for the sense recognition study (sec. IIVB[) and for the 
same reasons, the Maximum Gap method is getting less 
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FIG. 8: The median upper limit cross section <r me( j obtained 
with the three methods as a function of the angular reso- 
lution. The study is done in the case of 10 expected back- 
ground events and one expected WIMP event. The detector 
configuration is the following : sense recognition and energy 
threshold of 5 keV. 
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FIG. 9: Exclusion limits, WIMP-nucleon cross section as a 
function of the WIMP mass with the three methods and for 
three values of the energy threshold. The study is done in 
the case of 10 expected background events and one expected 
WIMP event. The upper bound of the energy range Er 2 is 
chosen as 200 keV. The detector configuration is the following 
: perfect angular resolution and sense recognition. 



competitive than the Poisson method beyond ere ~ 30°. 
Indeed, the expected WIMP distribution and background 
are getting closer and then, as previously explained, the 
Poisson method is more competitive than the Maximum 
Gap one. 



D. Effect of energy threshold 

As for direction-insensitive direct detection, the energy 
threshold plays a key role for directional detection. It 
depends on the target and the quenching factor. It is 
defined as the minimal recoil energy for which both the 
energy and the 3D track are well reconstructed. Reducing 
the energy threshold leads to a higher expected number 
of WIMP events and hence a more restrictive exclusion 
limit. Figure [9] presents the exclusion limits obtained, 
evaluated by the three statistical methods and for three 
different values of the energy threshold Er x , in the case of 
10 expected background events and one expected WIMP 
event. As the energy threshold Epj^ is taken between 5 
and 50 keV, the upper bound of the energy range Er 2 is 
chosen as 200 keV for this study. It can be noticed that 
the Likelihood method gives the best limit on the whole 
mass range and for the three threshold values. It is worth 
emphasizing that the effect of the energy threshold is very 
important even in the case of directional detection. Go- 
ing from a 5 keV to a 50 keV energy threshold leads to 
a loss of about one order of magnitude in exclusion lim- 
its. The situation is even worse at low WIMP masses. 
This outlines the fact that, as far as exclusion limits are 
concerned, lowering the energy threshold remains the ma- 
jor experimental issue for upcoming directional detector 
projects. 



E. Impact on limits from astrophysical 
uncertainties 

The effect of astrophysical uncertainties on exclusion 
limits obtained with direction-insensitive direct detection 
has been studied in details in [29|, [3(| ■ In the following, 
we investigate the effect on directional detection. 
Within the framework of a spherical isothermal halo 
model, two astrophysical parameters play a key role : 
the local dark matter density po and the local WIMP 
velocity dispersion vq. 

The first one is usually quoted within the range 
po ~ 0.2 - 0.8 GeV.c" 2 .cm- 3 and the value 0.3 GeV.c" 2 
is used as a " standard" value for the sake of comparison 
of various direct detector results HU . Recently, a value 
of the local Dark Matter density has been derived 
within the framework of a galaxy-model independent 
method (32[. The resulting local Dark Matter is 
po = 0.43 ± 0.11 ± 0.10 GeV/c 2 /cm 3 . As shown in QjJ, 
directional detection offers the possibility to constrain 
the WIMP-nucleon cross section, which can be relaxed 
into a constraint on po x a as a function of the WIMP 
mass. Then, the measured value of po [32| . together with 
its uncertainty, shall be accounted for when presenting 
discovery regions. However, when setting limits, the 
value of po does not change the shape of the exclusion 
limit, only its magnitude and with a negligible effect, 
owing to the orders of magnitude involved. In the 
following the so-called standard value is used. 

The second astrophysical parameter to be carefully 
handled is the WIMP velocity dispersion a v related 
to the asymptotic circular velocity vq, for which the 
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FIG. 10: Deviation of the median upper limit cross section 
(fmed) from the value a me d 220 obtained with the standard 
value vq = 220 km.s -1 , as a function of the WIMP mass. 
The study is done, for the three methods, with two extreme 
cases : 170 km.s -1 and 270 km.s - . 



"standard" value is 220 km.s -1 jU. As outlined in [29jj . 
recent determinations of its value span on a wide range 
and the impact for directionality has been studied, in 
the context of isotropy rejection (2?J ■ In the following, 
we investigate the effect of its value on directional 
exclusion limits by studying two extreme cases (170 and 
270 km.s- 1 ). 

Figure [TU] presents the deviation of the median up- 
per limit cross section (a me d) from the value er me£ ; 220 
obtained with the standard value vo = 220 km.s -1 
as a function of the WIMP mass. The study is done 
with a dataset of 10 background events and 1 WIMP 
event. For convenience, the Poisson case, represented 
by the black solid lines, is described first. As stated 
above, the Poisson exclusion limit is only sensitive to 
the total number of events and not to the directionality 
of the event angular distribution. When comparing 
low (170 km.s -1 , open circle) and high (270 km.s -1 , 
filled circle) values of vq, three effects on the excluded 
cross section can be noticed. The first one is linked 
to the energy threshold. For a low value of vq, the 
theoretical WIMP energy distribution is more peaked 
at low energy, leading to less events above the energy 
threshold. Hence, the cross section limit is higher in 
the case of vq = 170 km.s -1 than for vq — 270 km.s -1 . 
A second effect is due to the WIMP mass. Indeed, the 
difference between low and high values of vq is increased 
for low WIMP masses. Again, this is due a shift of a 
fraction of WIMP events below the threshold at low 
WIMP masses. 

When comparing with Maximum Gap and the likeli- 
hood methods, a third effect is introduced due to the 
directionality of the WIMP signal. The same general 
tendency is observed, but above a given WIMP mass 
(~ 40 GeV/c 2 ), the low vq exclusion limit becomes 



more restrictive. Indeed, a lower value of vq will induce 
a more directional WIMP event distribution as it is 
linked to the spread of the WIMP flux, see eq. [TJ Then, 
the difference between the background and the WIMP 
signal is enhanced, leading to a better exclusion. To 
summarize, when the WIMP mass is sufficiently high, 
the effect of the directionality at low vq overcomes the 
induced loss of events previously described. 
As a conclusion of this study, we can notice that the 
impact of uncertainties on the value of vq can lead to a 
deviation of the excluded cross section from 50% at low 
WIMP mass to 20% at large WIMP mass. 

As stated above, for this first study of an optimal exclu- 
sion method for directional detection, the standard halo 
model (spherical isothermal) has been considered. This is 
indeed the reference model to compare Dark Matter ex- 
periments and associated exclusion methods. Anisotropy, 
triaxiality, and clumping in the WIMP velocity distribu- 
tion would have an effect on directional detection sig- 
nal. As far as directional exclusion limits are concerned, 
this will mildly affect the result presented in this pa- 
per. Indeed, only the theoretical WIMP distribution used 
for the directional likelihood method would be changed, 
though remaining directional (anisotropic) and then dif- 
ferent from the background. However, the effect of non 
standard halo models (triaxial, with stochastic features 
or streams) will be addressed in a dedicated forthcoming 
paper [33j . 



V. PROSPECTING EXCLUSION LIMITS WITH 
DIRECTIONAL DETECTION 

To end-up this study, we present projected exclusion 
limits for a forthcoming directional detector proposed by 
the MIMAC collaboration [H]. We consider a 10 kg CF 4 
detector operated during ~ 3 years, allowing 3D track 
reconstruction, with a 10° angular resolution and a re- 
coil energy range 5-50 keV. In the following, a very low 
WIMP-nucleon cross section is considered, such as no 
WIMP event is expected during the acquisition time of 
the experiment. Note that a detector configuration has 
been chosen as an illustration, but the effect of the energy 
range, angular resolution has been shown previously. 
Figure Qj] presents the spin dependent cross section on 
proton (pb) as a function of the WIMP mass (GeV/c 2 ). 
Results are presented in the case of pure-proton approxi- 
mation [34| • The proton spin content has been chosen as 
(S p ) = 0.441 41] . Exclusion limits from direct detection 
searches are presented : COUPP S, KIMS [37j], NA- 
IAD UK, Picasso [HI andXenonlO Exclusion limits 
obtained with neutrino telescopes (Super K [Hj]) are also 
displayed. The theoretical region, obtained within the 
framework of the constrained minimal supersymmetric 
model, is taken from [35j . Constraints from collider data 
and relic abundance fl x h 2 , as measured with WMAP 5- 
year data [42], are accounted for. 
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FIG. 11: Spin dependent cross section on proton (pb) as a 
function of the WIMP mass (GeV/c 2 ). Results are presented 
in the case of pure-proton approximation (34|. The theoretical 
region, obtained within the framework of the constrained min- 
imal supersymmetric model, is taken from }35l ]. Constraints 
from collider data and relic abundance are accounted for. Ex- 
clusion limits from direct and indirect Dark Matter searches 
are displayed (see text for details). The projected exclusion 
limit of a forthcoming directional detector (30 kg.year) is pre- 
sented in three cases : background- free (solid line), with a 
background rate of 10 events/kg/year with sense recognition 
(dot-dashed line) and without sense recognition (dashed line) . 



Several cases are considered on figure [TT] : 

• the background-free measurement (dashed curve). 
It stands as the ultimate exclusion limit for a di- 
rectional detector with such an exposure. 

• with a background rate of 10 events/kg/year, which 
is large compared to current background measure- 
ments [3]. We consider the case of a directional 
detector with sense-recognition (dot-dashed line). 
The loss due to the presence of background is of 
the order of one order of magnitude, noticing that 
the use of the proposed Likelihood method allows 
the result to remain very satisfactory. 

• with a background rate of 10 events/kg/year but for 
a detector without sense- recognition (dashed line) . 



As stated above, as far as exclusion limits are con- 
cerned, the effect of not having sense-recognition is 
small. 

It highlights the exclusion power of a rather light direc- 
tional detector (10 kg), in the case that the discovery 
analysis [l(| failed at recognizing the galactic origin of 
the signal, for instance if the axial nucleon-WIMP cross 
section is very low. Even with highly contaminated data 
or without sense recognition directional detection with a 
10 kg CF4 detector would allow to set very constraining 
exclusion limits, about 3 orders of magnitude better than 
the one imposed by existing direct detectors. 



VI. CONCLUSION 

The conclusion of this study is threefold. 
First, a new likelihood method to derive exclusion limits 
has been proposed for the directional search of galactic 
WIMPs. Only the angular part of the spectrum is taken 
into account (in a given energy range), arguing that the 
energy part of the WIMP spectrum is featureless and may 
even be mimic by the background one. On the same real- 
istic set of directional data, the method has been shown 
to be more competitive than existing ones, or derived 
from existing ones. 

Second, in the spirit of directional detector opimization, 
a study of the effect of detector configurations on the 
exclusion limits has been done. Comparing the results 
from the three methods shows that the proposed ex- 
clusion Likelihood method gives the best results in all 
cases. Moreover, the effect of not having sense recogni- 
tion capability as well as having a low angular resolution 
has been shown to have little effect on exclusion lim- 
its. On the contrary, the effect of the energy threshold 
has been shown to largely influence the exclusion limits. 
This highlights the fact that the energy threshold is def- 
initely the major experimental challenge for upcoming 
directional detectors, while no sense recognition and low 
angular resolution could be sufficient as far as exclusion 
limits are concerned, which might be the case for low 
exposure data. Note that this conclusion does not hold 
for a discovery strategy 10], e.g. if the unknown axial 
WIMP-nucleon cross section lies in the 10~ 3 pb region. 
In this case, these detection issues remain of major inter- 
est, worth pursuing experimental efforts. 
Third, the exclusion power of a light directional detector 
(10 kg) has been shown to be competitive with respect 
to the existing exclusion limits. Even with highly con- 
taminated data or without sense recognition it allows to 
set exclusion limits about 3 orders of magnitude better 
than the one imposed by existing direct detectors. 
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